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Abstract. Non-Planckian (NP) spectral modifications of the CMB radiation spectrum can be pro¬ 
duced due to the existence of a non-zero value of the plasma frequency at the recombination epoch. 
We present here an analysis of NP effects on the cosmological radio background and we derive, for 
the first time, predictions of their amplitude on three different observables: the CMB spectrum, the 
Sunyaev-Zel’dovich (SZ) effect in cosmic structures, and the 21-cm background temperature bright¬ 
ness change. We find that NP effect can manifest in the CMB spectrum at v ^ 400 MHz as a drastic 
cut-off in the CMB intensity. Using the available CMB data in the relevant v range (i.e., mainly 
at ^ 1 GHz and in the COBE-FIRAS data frequency range), we derive upper limits on the plasma 
frequency v p = 206, 346 and 418 MHz at 1, 2 and 3 er confidence level, respectively. We find that 
the difference between the pure Planck spectrum and the one modified by NP effects is of the or¬ 
der of mJy/arcmin 2 at f ^ 0.5 GHz and it becomes smaller at higher frequencies where it is ~ 0.1 
mJy/arcmin 2 at v gj 150 GHz, thus indicating that the experimental route to probe NP effects in 
the early universe is to observe the cosmological radio background at very low frequencies. We have 
calculated for the first time the NP SZ effect (SZ np) using the upper limits on v p allowed by the 
CMB data. We found that the SZ np effect shows a unique spectral feature, i.e. a peak located 
exactly at the plasma frequency v p and this is independent of the cluster parameters (such as its 
temperature or optical depth). This offers therefore a way to measure directly and unambiguously 
the plasma frequency at the epoch of recombination by using galaxy clusters in the local universe, 
thus opening a unique window for the experimental exploration of plasma effects in the early universe. 
We have shown that the SKA-LOW has the potential to observe such a signal integrating over the 
central regions of high-temperature clusters. The studies of NP effects through the SZ np can be 
done by intensive observations of only one galaxy cluster, or with a stacked spectrum of a few well 
known clusters, thus avoiding the need of large statistical studies of source populations or wide area 
surveys. Finally, we also show that future low-u observations of the cosmological 21-cm background 
brightness temperature spectral changes have the possibility to set global constraints on NP effects by 
constraining the spectral variations of the temperature brightness change 5Tb induced by the plasma 
frequency value at the epoch of recombination. 
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1 Introduction 

The existence of a Planck spectrum for a system filled with a black body (BB) radiation assumes the 
presence of thermodynamic equilibrium between photons and matter. Interaction between photons 
and matter must be small enough for the gas of photons to be considered ideal and to avoid sub¬ 
stantial absorption and irreversible attenuation of electromagnetic radiation in the system (we note 
here that the interaction between photons is extremely weak). At the same time, the presence of 
matter and the weak interaction between matter and radiation is necessary for the photon gas to be 
in equilibrium [34] . These conditions are valid for a BB radiation from a photon gas and hence are also 
expected to be valid for the cosmic microwave background (CMB) spectrum emerging at the epoch 
of recombination (at redshift z rec ~ 1400) because the CMB spectrum is assumed and measured [22] 
to be close to a Planckian BB spectrum. 

In this context, the dispersion relation for electromagnetic radiation in non-ionized media (e.g., 
a rarefied gases) is given by 

oj = ck (1.1) 

where k is the wavenumber, c the speed of light and ui the angular frequency of the photon. However, 
in a medium which is ionized (e.g., an ionized plasma), the analogous dispersion relation is given by 
[48] 

to 2 = c 2 k 2 + n 2 p (1.2) 

where = 2nu p is the angular plasma frequency, with v p = Y2iLi riie 2 /Anmie being the actual 
plasma frequency, where the sum is over all charged particles species. The existence of such a disper¬ 
sion relation is due to the coupling between the electromagnetic waves (photons) and the collective 
behavior of the plasma [32]. Consequently, photons with frequencies lower than the plasma fre¬ 
quency v p do not exist in the final radiation spectrum. The presence of a non-zero value of the 
plasma frequency (in the following we use for our convenience the a-dimensional plasma frequency 
x p = hvp/ksTcMB, with ks being the Boltzmann constant and Tqmb the CMB temperature today) 
will cause the intensity of the photon spectrum to be zero for frequencies less than x p and this fact will 
reflect as a modified blackbody spectrum subject to non-Planckian (NP) effects. A general form for 
such a NP spectrum has been derived by [48] and [36] for astrophysical plasma-photon systems, such 
as the CMB and q-ray bursts atmospheres. We note here that the plasma frequency x p depends on 
the number density of all charged particles in the plasma, and hence deviation from Planck spectrum 
will also depend on this number density. 

In principle, the plasma frequency can be computed using knowledge of the cosmological recom¬ 
bination history. For example, a simple result obtained using a number density of 300 electrons/cm 3 
is that the plasma frequency at recombination is ss 0.2 MHz [48]. However, other effects bring un¬ 
certainties on this result. In fact, in addition to the ionized plasma, dark matter annihilation can 
produce electrons and positrons which will perturb the ionization fraction [6, 19, 25, 29, 39]. It has 
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been shown that the ionization fraction can increase due this effect by more than an order of magni¬ 
tude providing a plasma frequency of order of « 1.5 MHz [19]. Furthermore, electrons and positrons 
(hereafter we refer to these particles simply as electrons) can also be produced by hadronic collisions 
[17], originating another additional source of free particles, and consequently an increase of the plasma 
frequency value. 

The CMB spectrum is also vulnerable to other distortions like the /i and y distortions, related 
respectively to a frequency-dependent chemical potential and the Compton scattering of the CMB 
photons (see [47] for a review and references therein). The plasma effects in principle could also create 
a distortion well before the recombination era. Since the free particles number density depends on 
redshift, the plasma frequency for the CMB will also depend on redshift as v p oc (1 + z) 3 ^ 2 . Other 
distortions that can occur in the early universe are the Silk damping of small-scale perturbations, the 
cooling of photons by electrons and baryons, the decay and annihilation of relic particles, primordial 
magnetic fields and evaporating primordial black holes (see [7, 8, 47] and references therein). 

It is well known that some processes like Compton scattering, double Compton scattering and 
bremsstrahlung can change the photons energies or inject photons in the universe, in principle leading 
to a thermalization process, causing loss of information about the distortions occurring in the CMB 
[7, 42, 47]. However this does not seem to be possible for distortions due to the plasma frequencies. 
In fact, the time scales of these processes [47], t comp = 1.23 x 10 29 /(1 + z) 4 s, tdoubie = 1-34 x 
10 4o [a; 3 / (e x — 1)]/(1 + z ) 5 s and t/ ree = 8.59 x 10 26 [a; 3 / (e x — 1)]/(1 + z) 2 ' 5 s, are larger than the 
plasma time scale defined as t p = l/v p [38]. For example, at z = 1400 and for a plasma frequency of 
0.2 MHz corresponding to x p = 2.5 x 10 -9 , we obtain the values t comp = 3.2 x 10 16 s, tdoubie = 1.5 x 10 7 
s, tf ree = 72 s and t p = 6.4 x 10 -6 s; at z = 10 6 and a plasma frequency of 3 GHz corresponding 
to x p = 6.64 x 10~ 8 , we obtain t cornp = 1.2 x 10 5 s, tdoubie = 5.9 x 10 -5 s, tf ree = 3.8 x 10 -5 s and 
t p = 3.4 x 10 _1 ° s. The role of these processes, is, therefore, to eliminate the effects due to other 
distortions (see, e.g., [4, 5, 15, 16, 27, 28, 30, 31]), and to lead the CMB spectrum at the recombination 
to have the BB shape modified by the plasma frequency. This is the spectrum we can measure, and 
from which we can derive our results about the plasma frequency. 

So far, only theoretical arguments on NP plasma effects on the CMB have been discussed and 
neither experimental limit has been derived nor observational strategies to probe these effects has 
been proposed. All these uncertainties make it desirable to obtain limits to the plasma frequency 
from observations, that can be used to constrain the information about the processes occurred during 
and before the recombination. 

In this paper we want to derive for the first time constraints on the value of x p from the spectral 
shape of the CMB using current observations. In order to check if there is evidence for a non-zero 
value of Xp in the CMB spectrum, we therefore perform a comparison of the theoretical expectation 
for a NP CMB spectrum with the available observations at both low (radio) and high (microwave 
and mm.) frequencies. While the knowledge of the CMB spectrum is very precise in the microwave 
and mm. frequency regions after the COBE-FIRAS measurements [22], its knowledge at lower radio 
frequencies (particularly at frequencies below 1 GHz where the NP effects are expected to emerge) is 
still quite uncertain [26, 44, 45]. Therefore, in Sect.2 we will investigate first the limits on NP effects 
on the CMB spectrum and on the value of x p from the comparison with CMB measurements. 

We further stress that the modification on the CMB spectrum at \aw-v has an effect on other 
observable quantities that depend on the CMB spectrum. In this paper we will consider two additional 
cosmological observables whose spectral shapes are modified by the presence of NP effects in the CMB 
spectrum: i) the Sunyaev-Zel’dovich (SZ) effect produced by the inverse Compton scattering (ICS) of 
CMB photons off the electrons residing in the atmospheres of cosmic structures like galaxy clusters, 
radio galaxies and galactic halos [9], and ii) the temperature brightness change of the cosmological 21- 
crn background as produced by the evolution and reionization history of the primordial plasma during 
the Dark Ages (DA) and the epoch of reionization (EoR) [50]. Thus, in Sect.3 we study the imprint of 
the NP deviations on the CMB spectrum on the Sunyaev-Zel’dovich (SZ) effect in both thermal and 
non-thermal plasmas and we show that the main modifications to the standard SZ effect arise mainly 
in the very low-z/ range where a very specific spectral shape of the modified SZ effect can be observed. 
We will show that at v 500 MHz these spectral features allow to measure directly the plasma 
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frequency x p . Then in Sect.4 we will explore the effect on NP deviations of the CMB spectrum on the 
temperature brightness change of the cosmological 21-cm line background emerging from the Dark 
Ages and Epoch of Reionization and we predict the brightness temperature modifications observable 
with the next generation radio telescopes. We finally discuss our results in Sect.5. 

Throughout the paper, we use a flat, vacuum-dominated cosmological model with = 0.308, 
= 0.692 and H 0 = 67.8 km s -1 Mpc -1 [1]. 


2 The spectral distribution of the CMB 


The spectral distribution of a photon gas which is in equilibrium with matter (e.g rarefied gases) at 
temperature T is given by the Planck distribution which can be written in terms of the a-dimensional 
frequency x = hv/k B T as 


h(x) = 2 


ik B Tf 
(ihe ) 2 



( 2 . 1 ) 


For the case of the CMB, T — Tqmb = 2.725 ± 0.001 K [23] and h is the Planck constant, k B is the 
Boltzmann constant and c is the speed of light. In a ionized plasma with electron plasma frequency 
x p , deviations from Planck distribution in eq.(2.1) are possible because photons are suppressed at 
frequencies lower than x p . By taking in consideration this effect, we can write the generalized Planck 
distribution for the CMB photons as follows [48] 


Iq (x, x p ) — 2 


(■ k B T) 3 


(he) 


- 1 


\Jx 2 - x 2 H(ie - x p ), 


( 2 . 2 ) 


where H(x — x p ) is the Heaviside step function (with values 1 for x > x p and 0 for x < x p ). The 
step-function H(a; — x p ) takes into account the fact that photons are suppressed at frequencies below 
Xp. The difference between the two distributions is more pronounced if x p is larger than 1 but 
nonetheless there are appreciable differences even for values of x p < 1. The largest contrast between 
the generalized Planck distribution in eq.(2.2) and the standard Planck distribution in eq.(2.1) is 
that the former has a cut-off frequency which depends on x p . We also notice that the frequency at 
which the maximum of the CMB spectrum occurs shifts to higher frequencies as the cuff-frequency 
x p increases [48] and that the CMB peak intensity decreases as the value of x p increases. 

In order to investigate the presence of non-Planckian effect on the CMB spectrum, we fit the 
NP distribution of the CMB spectrum in eq.(2.2) to the CMB data [22],[26, 44, 45] and we minimize 
the x 2 with respect to the plasma frequency x p \ in principle one can also use the CMB temperature 
T as a free parameter: however, we set its value to the one measured by COBE-FIRAS that is very 
accurate with uncertainty of order of « 0.03% and it is also determined in the high -v range of the 
CMB spectrum (around its maximum) where the NP effects are small. With this fitting procedure 
we obtain a value of x p = 0 for the minimum y 2 . This means that a pure Planckian spectrum is best 
fitting the present data. However, the uncertainties in the data (in particular at low frequencies v < 1 
GHz) allow us to set upper limits on the quantity x p which are x p = (3.63, 6.10, 7.36) x 10~ 3 for lcr, 
2 a and 3er confidence level (c.l.), corresponding to frequencies of 206, 346 and 418 MHz, respectively. 
In Fig 1 we show the NP effects on the CMB spectrum over the full frequency range. This result 
indicates that it is not possible at the 3cr c.l. to distinguish between the Planck spectrum and the 
NP modified spectrum based on the current knowledge of the CMB spectrum. This fact opens hence 
an experimental window for the exploration of NP effects in the early universe. We then show in Fig 
2 the CMB spectrum and its experimental uncertainties zooming around the low frequency region at 
v < 1 GHz showing the modified CMB spectrum with the upper limits on x v at lcr, 2 a and 3cr c.l. 
These results show that the intensity of the CMB goes to zero below ss 206 MHz, i.e. for the value of 
x p = 3.63 x 10~ 3 at lcr c.l., or below « 418 MHz, i.e. for the value of x p = 7.36 x 10 -3 at 3cr c.l. In 
Fig 3, we show the relative CMB spectrum variations around the peak of the CMB where the COBE 
data [22] are available (for the sake of illustration of the amplitude of the NP effects we zoom within 
the region 149-191 GHz around the maximum of the CMB spectrum). As shown in Fig 3, the CMB 
spectrum intensity decreases with increasing values of x p and the frequency location of its maximum 
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Figure 1 . The non-Planckian spectral distribution of the CMB for different values of the plasma frequency 
u p as derived from the fit to the data. The experimental data shown here are from [26, 44, 45] and COBE [22]. 
Other data in the range ~ 1.3 — 50 GHz obtained from ground-based, ballon-borne and from the COBE-DR 
experiment (see http://asd.gsfc.nasa.gov/archive/arcade/cmb_spectrum.html) are not shown here but they lie 
almost exactly along the curves of the CMB spectrum shown in the plot. 


increases with increasing value of x p . The amplitude of the CMB intensity variation due to NP effects 
is of order of ss mJy/arcmin 2 at low-i/ while it is of order of ss 0.1 mJy/arcmin 2 at high-i/. Therefore, 
the best frequency region to look experimentally for these effects is at low radio frequencies, i.e. at 
v 1000 MHz. Our analysis shows that precise measurements of the CMB spectrum at low radio 
frequency have the power to set strong constraints to NP effects due to the existence of a finite value of 
the plasma frequency at recombination and that one can probe these effects with the next generation 
experiments like the SKA (SKA-LOW and SKA-MID [18]) and HERA [2], 

Even though the CMB spectrum can be considered as a direct probe of NP effects in the early universe, 
its measurement at very low radio frequencies 400 MHz is challenging and it is subject to various 
experimental systematics and biases, and on the ability to subtract carefully radio foregrounds [37]. 
The next generation low-r' radio telescopes (SKA, HERA) will increase our knowledge of the CMB 
spectrum in this frequency range, but the problem of foreground contamination and of the component 
separation of such diffuse signals at low frequency still remains. Therefore it is interesting to explore 
the use of other observable quantities that depend directly on the CMB spectrum, like, e.g., the 
modified SZ effect [11] and the modifications to the 21-cm brightness temperature change w.r.t. the 
CMB. In the following Sections we will discuss these two additional probes of NP effects in the early 
universe. 
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Figure 2. The non-Planckian effects due to a finite value of the plasma frequency v p on the CMB spectrum 
at v < 1 GHz for different values of the upper limit on v v as obtained from the fit to the CMB spectrum data. 
Experimental data are from [26, 44, 45]. 


3 The modified Sunyaev ZePdovich effect 

Another method which can be used to probe the existence of non-Planckian effect in the CMB is to 
use the modified SZ effect (see [11], and [3, 10] for the general derivation of the SZ effect). The SZ 
effect is a spectral distortion of the CMB spectrum where the CMB photons interact with the electrons 
(thermal and non-thermal) residing in the atmosphere of galaxy clusters and other structures on large 
scales, like e.g. radio galaxies and galaxy halos. The SZ effect is a differential measure of the CMB 
spectrum on and off an area of the sky containing the selected cosmic structure and hence it is not 
affected by large-scale foregrounds in the observations at low-frequency (see [13] for a discussion). 

The spectral distortion due to SZ effect on the CMB spectrum can be computed, at first order 
in the optical depth r e = <jt f d£n e , as follows [10]: 


A I(x) = T e 


dsP 1 (s)I 0 {xe s ) — Iq(x) , 


(3.1) 


where -Pi(s) is the single scattering redistribution function, obtained by summing the photon redis¬ 
tribution function P s (s,p), that gives the probability to have a logarithmic shift s = ln(z//i/) in the 
photon frequency due to the ICS process by an electron with a-dimensional momentum p = /3y, 
over the electron momentum distribution function / e (p), normalized as to have f^° f e (p)dp = 1 (see 
[10, 20]). This is written as 

/»00 

P 1 (s)= dpf e (jp)P s {s,p). (3.2) 

Jo 

For the calculation of the standard SZ effect, the quantity Iq(x) is a Planck spectrum in the case of 
the CMB. Since the SZ effect is related directly to the spectrum of the CMB (see eq.3.1), the existence 
of non-Planckian effect in the CMB background will be reflected directly in the SZ spectral distortion, 
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Figure 3. The effect of non-Planckian distribution on the CMB spectrum in frequency range 143 GHz to 
191 GHz is shown for different values of the plasma frequency x p as derived in Sect.2. Experimental data are 
from COBE [22], 


and it is expected to be visible in particular around the cut-off frequency x p . In order to compute the 
SZ effect for the case of a NP spectrum, eq.(2.2) is inserted into eq.(3.1) instead of 7 0 (x) as follows 


A I{x] x p ) = 


dsPi(s)Io(xe s -,x p ) - I 0 (x;x p ) 


(3.3) 
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Figure 4. The thermal SZ effect spectral distortions computed for galaxy clusters with increasing plasma 
temperature (see various panels) for a standard Planck distribution (solid line) and including the effect of 
non-Planckian distribution of photons for the values of the plasma frequency x v derived at la, 2a and 30- 
level. The cluster plasma optical depth is fixed to the value r = 0.001. 


One can notice directly from this equation that the final SZ spectral distortion in the case of NP 
plasma effect depends directly on x pi hence on the plasma frequency u p . 

In the context of the SZ effect modified by NP effects (that we refer here to as SZ^rp) we first 
compute the thermal spectral distortion for an electron plasma within galaxy clusters with r e = 10~ 3 
and with different electron temperatures of 5 keV, 10 keV, 15 keV and 20 keV. We show in Fig 4 
the SZjvp spectrum for different values of x p derived from the analysis of the CMB spectrum in the 
previous Sect.2. 

We notice that the frequency at which the maxima of the SZ^vp occurs is exactly x p and is 
independent of the cluster temperature. This spectral feature offers therefore a direct experimental 
probe of the plasma frequency in the early universe as measured by the SZ jvp in structures of the local 
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Figure 5. The non-thermal SZ spectral distortion for the case of non-thermal plasmas with increasing electron 
minimum momentum p\ (see various panels) for a usual Planck distribution (solid line) of photons and the 
effect of non-Planckian distribution at lcr, 2a and 3<r level computed in the case of a single power law for 
different minimum momentum p\. We use here an optical depth of r = 0.001. 


universe. We stress that this is a unique characteristic of the SZ mp and allows to use observations of 
the local universe to infer directly the physics of the early universe. 

We also show that the amplitude of the peak of SZjyp increases with the temperature and density 
(i.e. with the total Compton parameter Y oc f d£P ei where P e is the total electron pressure) of the 
electron gas. The spectral shape of the thermal SZjyp is very peculiar and it is peaked in a quite 
narrow frequency range of the order of Au a; 20 MHz for kT = 5 keV and Av « 50 MHz for kT = 20 
keV, reflecting hence the relativistic effects of the electron distribution in the photon re-distribution 
function P\{s). At frequencies larger than 1 GHz, the SZjyp becomes close to the standard Planckian 
one, i.e. the standard SZ effect calculated with the non-distorted CMB Planck spectrum [10] as an 
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Figure 6. The difference between the SZjvp and the Planckian SZ for a Bullet-like cluster with temperature 
of 15 keV (see [35]) subtending an angle of 5' to its Rsoo- An optical depth of r = 0.001 is assumed. The 
SKA1 (red) and SKA-50% (orange) sensitivity is calculated for 1000 hrs integration while the eVLA sensitivity 
(green) is calculated for 12 hrs integration. 


SKA-50% 


\ i\ \ / ; 

SKA-1 • 


.v p =206 MHz 

. Vp=346 MHz 
i/ p =418 MHz 


eVLA 


• ••IM MIAIAIM MUIXIUMIAIA 


input. 

The SZjvp can appear in the form of both thermal and non-thermal effects depending on the 
nature of the electron population producing the ICS of the CMB photons (i.e. thermal electrons in 
clusters for the thermal SZjvp and non-thermal or relativistic electrons for the non-thermal SZjvp). 
For completeness, therefore, we also compute the effect of non-planckian CMB distribution on the 
non-thermal SZjvp effect which is shown in Fig 5 for a power law electron spectrum with spectral index 
a = 2.5 and different values of the minimum momentum pi of the non-thermal electron distribution. 
In this case, the peak of the SZjvp decreases with increasing value of pi because the high-E electrons 
scatter photons to very high frequency. In this case the peak of the SZjvp is less pronounced (see 
Fig. 5) but it is still located at the plasma frequency value x p . The width of the SZjvp spectrum (see 
Fig. 5) is larger due to the enhanced impact of the relativistic effects of the high-E electron population 
and its shape reflects therefore the different nature of the scattering plasma. Possible observations of 
the SZjvp can therefore address the question of the intrinsic nature of the plasma in the target cosmic 
structure (e.g., galaxy clusters vs. radio galaxies). 

In order to study the possible detectability of NP effect on the CMB, we have compared the SZjvp 
spectrum integrated over an angle of 5 7 radius of a Bullet-like cluster with temperature of kT =15 
keV. Fig 6 shows the difference between the SZ spectrum A I and A / compared to the sensitivity 
of the SKA1 and SKA50% for 1000 hrs of observations and to the eVLA sensitivity for 12 hrs of 
observation. We found that the values on the plasma frequency v p = 206, 346 and 418 MHz at 1, 2 
and 3 a c.I, derived from the analysis of the CMB spectrum, can be detected with both SKA-LOW 
and SKA-MID bandl (350-1050 MHz). 
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4 Modifications to the cosmological 21-cm background 


The frequency range where NP effects can be observed (see Sects.2 and 3) is the same where the 
distortion of the CMB spectrum induced by the 21-cm background is more relevant [24]. For this 
reason we explore in this Section how NP effects on the CMB spectrum as in eq.(2.2) impact on the 
cosmological 21-cm background spectrum. To this aim, we discuss the results of our spectral analysis 
in term of the brightness temperature change given by 


5T{v) 


2k bv 2 


81{v). 


(4.1) 


For a Planckian CMB spectrum, the 21-cm background appears as a perturbation on the CMB 
background, so we have: 

5T(u)=T 21 (u)-T 0 (u), (4.2) 

where X 21 is the brightness temperature of the 21-cm background and To is the brightness temperature 
of the CMB. In the case of a non-planckian CMB spectrum, the previous relation can be generalized 
as follows 

5T(v\ v p ) = T 2 i(zq i/ p ) - T 0 (t/; v p ). (4.3) 

Once we fixed the value of x p , we can distinguish two frequencies regimes: 

i) v < v p : in this case T 0 = 0, and therefore ST = T 2 1- If we assume that T 2 1 « T 2 i, i.e. the 21-cm 
background is not sensibly changed by the NP spectral distortion, as it is expected because the main 
physical processes affecting the 21-cm background depend on the global temperature of the system 
which is not heavily affected by a distortion at small frequencies, then we can write: 


ST = ST + T 0 , 


(4.4) 


so the frequency change of the brightness temperature is the same as in the Planckian case, but its 
amplitude is shifted by a value T 0 (that in the RJ region has a constant value equal to the CMB 
temperature); 

ii) v > v v \ in this case we can write the NP modified spectrum in the RJ region as: 


To oc X\/x 2 — Xp 


and the corresponding brightness temperature is: 




By assuming again T 2 i « T 2 i, the change in the temperature brightness is given by: 

X 2 j ' 


6T = ST + T 0 -T 0 = 6T + T 0 1 -\ 1 


(4.5) 


(4.6) 


(4.7) 


In Fig. 7 we show the spectral change of ST for different values of v p from 0.2 to 100 MHz, by assuming 
a standard benchmark model for ST (as discussed in[12, 14]). For v < v v the frequency change of ST 
is the same ST obtained in the standard Planckian case but its amplitude is increased by the value 
To; for v 33 v p the quantity ST decreases very rapidly with increasing frequency, and for v v p it 
then becomes equal to the standard value ST in the Planckian case. Fig. 7 shows that even low values 
of the plasma frequency v p 33 206 MHz (i.e., the 1 <7 upper limit set by the CMB spectrum analysis 
in Sect.2) can largely change the 21-cm radiation background spectrum and hence our ability to 
recover the history of the EoR and of the DA. NP effects do not change the physical mechanisms that 
are working during these epochs but change the brightness temperature contrast ST relative to the 
CMB that we can use to prove these remote cosmic epochs. Interestingly, the brightness temperature 
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Figure 7. The change in brigthness temperature ST for v p = 0 (solid line), v p = 0.2 (dotted), 1 (dashed), 10 
(dot-dashed), 20 (three dots-dashed) and 100 MHz (long dashed) as a function of the frequency. 


contrast ST in the case in which NP effects are present on the CMB is increased w.r.t. the standard 
assumption of a pure Planckian CMB spectrum. This fact on one hand yields larger signals that could 
be detected more easily even with precursors of the SKA and HERA, and on the other hand provides 
an effective possibility to set limits or even prove the existence of NP effects in the early universe by 
looking at the global spectrum of the 21-cm background. Available limits from the PAPER experiment 
[40] in the 100-200 MHz provide limits only on the average temperature brightness of (T&) < 275 and 
291 mK for values of the ionization fraction Xi = 0.5 and Xi = 0.75 of the ionization power spectra, 
respectively, which cannot be unfortunately directly compared to the predictions for the quantity ST 
used in our paper. The next generation SKA and HERA low-frequency interferometers will have the 
possibility to improve the knowledge on the brightness temperature contrast ST down to levels of 
~ 3 mK (~ 0.3 mK) with SKA phase-1 (SKA phase-2) at 3cr c.l. in about 1000 hrs at 150 MHz [33] 
and possibly cover a wider frequency range, thus allowing to set much stronger limits on v p from the 
spectrum of ST(v). 

The SKA1-LOW band 50-350 MHz will be able to set limits only on plasma frequency values 
v p 10 MHz because for lower values of v p of order of a few MHz the NP modified 21-cm background 
spectrum is very close to the Planckian case at the lower end of the SKA-LOW band (see Fig. 7). 
Therefore, the SKA-LOW will have a sensitivity to the NP effects at the redshift (frequency) of the 
EoR. We speculate here that other sources of more recent NP effects can arise during the EoR due 
to the increasing ionization fraction of the cosmic plasma at this epoch caused by Ly-a interactions, 
heating due to UV/X-ray primordial sources and DM annihilation. The treatment of these secondary 
NP effects is well beyond the purposes of this initial paper and will be addressed more specifically 
elsewhere. 
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5 Discussion and conclusions 


We have performed for the first time in this paper an analysis of NP effects on the cosmological radio 
background and we derived predictions for their amplitude on three different observables: the CMB 
spectrum, the SZ effect and the 21-cm temperature brightness change. 

We have shown that NP effect can manifest on the CMB spectrum at v ;$ 400 MHz as a drastic 
cut-off in the intensity of the CMB (see Fig. 2). Using the CMB data of COBE and lower frequency 
measurements, we derived for the first time upper limits on the plasma frequency v p = 206, 346 and 
418 MHz at 1, 2, 3 a c.l., respectively. The difference between the pure BB Planck spectrum and the 
one modified by NP effects at these low frequencies is of the order of mJy/arcmin 2 and this difference 
becomes smaller at higher frequencies (v « 150 GHz) where it is of the order of 0.1 mJy/arcmin 2 , 
thus indicating that the experimental route to probe NP effects in the early universe is to observe the 
cosmological radio background at very low frequencies. 

We have also computed the SZ^p effect using the upper limits on v p allowed by the present 
data of the CMB (Sect.2) and it has been found that the NP effect produces a unique imprint on the 
SZjvp effect spectrum at v 1 GHz, i.e., a peak located exactly at the plasma frequency v p and that 
this is independent of the cluster parameters (such as its temperature or optical depth). This offers a 
way to measure directly and unambiguously the plasma frequency in the early universe at the epoch 
of recombination by using (in principle only one) cluster in the local universe, thus opening a unique 
window for the experimental exploration of plasma effects in the early universe. We have shown that 
the SKA-LOW and SKA-MID has the potential to observe such a signal integrating over the central 
regions (~ 5' radius) of high-temperature ( k B T ~ 15 keV) clusters. 

Observing the SZ np also benefits from its differential nature thus being less affected by the large 
impact of large-scale foreground emission that is one of the main systematic biases that limits the 
study of the intensity spectrum of the cosmological background radiation. The SZjyp has also the 
appealing property that we can study the presence of NP effects in the early universe by looking 
at very local cosmic structures for which the structural parameters are known with high accuracy. 
Finally, we mention that studies of NP effects through the SZ^rp can be done by intensive observations 
of only one galaxy cluster, or with a stacked spectrum of a few well known clusters, thus avoiding the 
need of large statistical studies of source populations or wide area surveys. 

Finally, we have shown that future low-frequency observations of the cosmological 21-cm bright¬ 
ness temperature spectral changes have the possibility to set global constraints on NP effects by 
constraining the spectral variations of ST induced by the plasma frequency value at the epoch of 
recombination. We discussed, in this context, that even moderate limits on the average brightness 
temperature of the 21-cm background obtained with SKA precursors, like e.g. the PAPER experi¬ 
ment, are able to start limiting the possible values of v p in its high-frequency domain (of order of 100 
to a few hundreds MHz). 

In conclusion, we have demonstrated that the study of low-frequency cosmological radio back¬ 
ground has a strong and unique potential for proving the physics of the early universe. This paper 
aligns such the idea of proving NP effects with our previous studies of the photon decay effects in the 
early universe [11] and on the study of the DA and EoR through the SZE-21cm [12, 13], and indicates 
that this area of investigation of the fundamental physics of the universe will receive a boost with the 
next generation high-sensitivity radio telescopes like the SKA. 

Acknowledgments 

S.C. acknowledges support by the South African Research Chairs Initiative of the Department of 
Science and Technology and National Research Foundation and by the Square Kilometre Array (SKA). 
P.M. and M.S.E. acknowledge support from the DST/NRF SKA post-graduate bursary initiative. 


References 

[1] Ade, P.A.R., et al., 2015, arXiv: 1502.01589 


- 12 - 


[2] Backer, D.C.et al. 2010, Submitted for consideration by the Astro2010 Decadal Survey Program Panel 
RMS: Radio (Meter/Centimeter) and Millimeter/Submillimeter 

[3] Birkinshaw, M. 1999, Phys. Rep., 310, 97 

[4] Burigana Cl, Danese L., de Zotti G., 1991a, ApJ, 379, 1 

[5] Burigana C., Danese L., de Zotti G., 1991b, A&A, 246, 49 

[6] Chen, X. & Kamionkowski, M., 2004, Phys.Rev. D, 70, 043502 

[7] Chluba, J. & Sunyaev, R.A. 2012, MNRAS, 419, 1294-1314 

[8] Chluba, J. 2014, arXiv:1405.6938vl [astrop-ph.CO] 

[9] Colafrancesco, S. 2013, Acta Pol., 53, 560 

[10] Colafrancesco, S., Marchegiani, P. and Palladino, E., 2003, A&A, 397, 27-52 

[11] Colafrancesco, S., and Marchegiani, P., 2014a, A&A, 562, L2 

[12] Colafrancesco, S., and Marchegiani, P., 2014b, A&A, submitted 

[13] S. Colafrancesco, P. Marchegiani and M.S. Emritte: Cosmology with galaxy clusters: studying the Dark 
Ages and the Epoch of Rcionization in the SKA era, 2015, in proceedings of ’’Advancing Astrophysics 
with the Square Kilometre Array”, PoS(AASKA14)034, in Press 

[14] Cooray, A., 2006, Phys. Rev. D, 73, 103001 

[15] Danese L., de Zotti G., 1977, Nuovo Cimento Riv. Ser., 7, 277 

[16] Danese L., de Zotti G., 1982, A&A, 107, 39 

[17] Dernier, C.D., 1986, ApJ, 307, 47 

[18] Dewdney, P., Turner, W., Millenaar, R., McCool, R., Lazio, J. & Cornwell, T., 2012, SKA baseline 
design document, http://www.skatelescope.org/wp-content/uploads/2012/07/SKA-TEL-SKO-DD-001- 
1 -BaselineDesignl.pdf 

[19] Dvorkin, C., Blum, Kfir. & Zaldarriaga, M. 2013, Phys. Rev. D, 87, 103522 

[20] Enfilin, T.A. & Kaiser, C.R., 2000, A&A, 360, 417 

[21] Evoli, C., Mesinger, A. & Ferrara, A., 2014, arXiv: 1408.1109 

[22] Fixsen, D.J. et al., 1996, ApJ, 473, 576 

[23] Fixsen, D.J. & Mather, J.C. 2002, ApJ, 581, 817 

[24] Furlanetto, S.R., Oh, S.P., & Briggs, F.H., 2006, Phys. Rep., 433, 181 

[25] Galli, S. et al., 2009, Pliys.Rev.D, 80, 023505 

[26] Howell, T.F. & Shakeshaft, J.R. , 1967, Nature, 216, 753 

[27] Hu W., Silk J., 1993a, Phys. Rev. D, 48, 485 

[28] Hu W., Silk J., 1993b, Phys. Rev. Lett., 70, 2661 

[29] Hiitsi, G., Hektor, A. & Raidal, M., 2009, A&A, 505, 999-1005 

[30] Illarionov A. F., Sunyaev R. A., 1975a, SvA, 18, 413 

[31] Illarionov A. F., Sunyaev R. A., 1975b, SvA, 18, 691 

[32] Kittell, C. 1986, Introduction to Solid State Physics 

[33] Kopmans, L. 2010, 

http://www.mrao.cam.ac.uk/projects/aa.vp/presentations/Koopmans-EoRSKAl.pdf 

[34] Landau, L.D. and Lifshitz, E. M. Statistical Physics 

[35] Markevitch, M., G. et al. 2002 ApJ, 567, L27 (Nauka, Moscow, 1976; Pergamon, Oxford, 1980). 

[36] Medvedev, M.V. 1999, Phys.Rev.E, Vol 59, 5 

[37] Oliveira-Costa, A. et al. 2008, MNRAS, 388, 247-260 


- 13 - 



[38] Padmanabhan, N. 2000, ’’Theoretical Astrophysics”, Vol 1, Cambridge University Press 

[39] Padmanabhan, N. & Finkbeiner, D., P., 2005, Phys.Rev.D, 72,023508 

[40] Parsons, A.R. et al., 2014, ApJ, 788, 106 

[41] Planck Collaboration 2011, A&A, 536, 28 

[42] Procopio, P. & Burigana C., 2009, A&A, 507, 1243 

[43] Rogers, A.E.E. & Bowman, J.D. 2008, ApJ, 136, 641-648 

[44] Sironi, G. et al. 1990, ApJ, 357, 301 

[45] Sironi, G. et al. 1991, ApJ, 378, 550 

[46] Smoot, G., F 1997 astro-pli/9705101v2 

[47] Tashiro, H. 2014, Prog.Theo. Exp. Phys, 06B107, DOI: 10.1093/ptep/ptu066 

[48] Triger, S.A. & Khomkin, A.L., 2010, Plasma Physics Reports, 36, 1095 

[49] Valdes, M. et al., 2013, MNRAS, 429, 1705 

[50] Zaroubi, S. 2013, in The First Galaxies, Astrophysics and Space Science Library, 396, 45 


- 14 - 



